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In [1] : # Defineerime impordid, konstandid ja 

# abimeetodid edaspidiseks kasutamiseks 

import matplotlib.pyplot as pit 
import numpy as np 

pi = np.pi 
sqrt = np.sqrt 
e = np . e 

def fig (x=16, y=14) : 

pit . figure ( figsize= (x, y) ) 

pit . rcParams [' font . size' ] = 18 

def get_x(a, pm, step=0.001) : 

return np . arange (a-pm, a+pm, step) 



Ulesanne 71. a) 



Leidke funktsiooni f(x) = ^/x Taylory n-jarku poliinoom punktis a, kus a = 1, n = 2, 



In [2] : a = 1 

def yl71a (x) : 

return sqrt (x) 

def yl71a2 (x) : 

return 1 + (x - a) 12 - ( (x - a)**2)/8 

xl = get_x (a, 1 ) 
x2 = get_x ( 5 , 5 ) 

fig() 

pit . figure ( 1 ) 
pit . subplot (221) 

plt.plot(xl, yl71a(xl), ' r' , label=' $ \ sqrt { x } $ ' ) 
plt.plot(xl, yl71a2(xl), ' g' , label=' $P (x ) $ ' ) 
pit. plot (a, yl71a(a), 'ro') 
pit . legend (loc=' best' ) 

pit . subplot (222) 

pit. plot (x2, yl71a(x2), ' r' , label=' $ \ sqrt { x } $ ' ) 
pit. plot (x2, yl71a2(x2), ' g' , label = ' $P (x ) $ ' ) 
pit. plot (a, yl71a(a), 'ro') 
pit . legend (loc=' best' ) 



pit . show ( ) 
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Ulesanne 71. b) 



Leidke funktsiooni f(x) = Taylory rt-jarku poliinoom punktis a, kus a = 25, n = 2. 



In [3] : yl71b = yl71a 

a = 25 



def yl71b2 (x) : 

return 5 + (x - a)/10 - ( (x - a)**2)/1000 

xl = get_x (a, 10) 
x2 = get_x(a+25, 50) 

fig() 



pit . figure ( 1 ) 
pit . subplot (221) 

plt.plot(xl, yl71b(xl), 'r', label=' $ \ sgrt { x } $ ' ) 
pit. plot (xl, yl71b2(xl), 'g', labels' $P (x) $' ) 
pit. plot (a, yl71b(a), 'ro') 
pit . legend (loc=' best' ) 

pit . subplot (222) 

pit. plot (x2, yl71b(x2), ' r' , label=' $ \ sgrt { x } $ ' ) 
pit. plot (x2, yl71b2(x2), ' g' , label = ' $P (x) $ ' ) 
pit. plot (a, yl71b(a), 'ro') 
pit . legend (loc=' best' ) 



pit . show ( ) 




Ulesanne 71. c) 

Leidke funktsiooni f(x) = arctanx Taylory n-jarku poliinoom punktis a, kus a = 1, n = 3. 



In [4] : a = 1 

def yl71c (x) : 

return np.arctan(x) 

def yl71c2 (x) : 

return pi/4 + (x - a) /2 - (x - a)**2/4 + (x - a) **3/12 

xl = get_x (a, 1 ) 
x2 = get_x (a, 3 ) 

fig() 

pit . figure ( 1 ) 
pit . subplot (221) 

plt.plot(xl, yl71c(xl), ' r' , label=r ' $ \arctan { (x) } $ ' ) 
plt.plot(xl, yl71c2(xl), 'g', label = ' $P (x ) $ ' ) 
pit. plot (a, yl71c(a), 'ro') 
pit . legend (loc=' best' ) 

pit . subplot (222) 

pit. plot (x2, yl71c(x2), ' r' , label=r ' $ \arctan { (x) } $ ' ) 
pit. plot (x2, yl71c2(x2), ' g' , label = ' $P (x ) $ ' ) 
pit. plot (a, yl71c(a), 'ro') 
pit . legend (loc=' best' ) 



pit . show ( ) 
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Ulesanne 71. d) 



Leidke funktsiooni f(x) — s'mx Taylory n-jarku poliinoom punktis a, kus a = ^, n = 4. 

In [5] : a = pi/6 



def yl71d (x) : 

return np.sin(x) 



def yl71d2 (x) : 

d = lambda x : x - a 

return 1/2 + (sqrt (3) /2) *d(x) \ 

- (1/4) *d(x) **2 - (sqrt (3) /12) *d(x) **3 \ 
+ (1/48) *d (x) **4 

xl = get_x (a, pi) 

x2 = get_x(a, 1.5*pi) 

fig() 



pit . figure ( 1 ) 
pit . subplot (221) 

plt.plot(xl, yl71d(xl), 'r', label = ' $ \ sin { (x ) } $ ' ) 
pit. plot (xl, yl71d2(xl), 'g', label=' $P (x) $ ' ) 
pit. plot (a, yl71d(a), 'ro') 
pit . legend (loc=' best' ) 

pit . subplot (222) 

pit. plot (x2, yl71d(x2), ' r' , label=' $ \ sin { (x ) } $ ' ) 
pit. plot (x2, yl71d2(x2), ' g' , label = ' $P (x) $ ' ) 
pit. plot (a, yl71d(a), 'ro') 
pit . legend (loc=' best' ) 



pit . legend (loc=' best' ) 



pit . show ( ) 
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Ulesanne 73. a) 



Hinnake viga Taylori valemi kaudu saadud ligikaudses valemis: sinx sa x — kui \x\ < |. 

In [6] : def yl73a (x) : 

return np.sin(x) 

def yl73a2 (x) : 

return x - x**3/6 

xl = get_x(0, 0.5) 
x2 = get_x(0, 4) 

fig() 

pit . figure ( 1 ) 
pit . subplot (221) 

plt.plot(xl, yl73a(xl), ' r' , label=' $ \ sin { (x ) } $ ' ) 
plt.plot(xl, yl73a2(xl), ' g' , label=r' $x - \ f rac { x"3 } { 6 } $ ' ) 
pit . legend (loc=' best' ) 

pit . subplot (222) 

pit. plot (x2, yl73a(x2), ' r' , label=' $ \ sin { (x ) } $ ' ) 

pit. plot (x2, yl73a2(x2), ' g' , label=r' $x - \f rac{x A 3 } { 6 } $' ) 

pit . axvline (x=-0 . 5, color='k', ls=' dashed' ) 

pit . axvline (x=0 . 5, color=' k' , ls='dashed') 

pit . legend (loc=' best' ) 



pit . show ( ) 
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Ulesanne 73. b) 

Hinnake viga Taylori valemi kaudu saadud ligikaudses valemis: tana; 

In [7] : def yl73b (x) : 

return np.tan(x) 

def yl73b2 (x) : 

return x + x**3/3 

xl = get_x(0, 0.1, step=0.001) 
x2 = get_x(0, 1.5, step=0.001) 



fig() 

pit . figure ( 1 ) 
pit . subplot (221) 
plt.plot(xl, yl73b(xl), 
pit. plot (xl, yl73b2(xl) 
pit . legend (loc=' best' ) 



kui Id < 0.1. 



label=r' $\tan{ (x) }$' ) 
label=r'$x + \f rac{x A 3 } { 3 } $' ) 



pit . subplot (222) 

pit. plot (x2, yl73b(x2), ' r' , label=r ' $ \tan { (x) } $ ' ) 

pit. plot (x2, yl73b2(x2), 'g', label=r' $x + \f rac{x A 3 } { 3 } $' ) 

pit . axvline (x=-0 . 1, color='k', ls=' dashed' ) 

pit . axvline (x=0 . 1 , color=' k' , ls='dashed') 

pit . legend (loc=' best' ) 



pit . show ( ) 




Ulesanne 73. c) 



i 

Hinnake viga Taylori valemi kaudu saadud ligikaudses valemis: yl + a; ~ 1 + | — kui 0 < a; < 1. 



In [8] : def yl73c (x) : 

return sqrt(l+x) 

def yl73c2 (x) : 

return 1 + x/2 - (x**2)/8 

xl = get_x(0.5, 0.5) 
x2 = get_x ( 2 , 2 ) 

fig() 

pit . figure ( 1 ) 
pit . subplot (221) 

pit. plot (xl, yl73c(xl), 'r', label=' $ \ sqrt { 1+x } $ ' ) 

plt.plot(xl, yl73c2(xl), ' g' , label=r ' $ l+\ f rac { x } { 2 } -\ f rac { x A 2 } { 8 } $ ' ) 
pit . legend (loc=' best' ) 

pit . subplot (222) 

pit. plot (x2, yl73c(x2), ' r' , label=' $ \ sqrt { 1+x }$' ) 

pit. plot (x2, yl73c2(x2), ' g' , label=r ' $ l+\ f rac { x } { 2 } -\ f rac { x A 2 } { 8 } $ ' ) 
pit . axvline (x=0 , color='k', ls=' dashed' ) 
pit . axvline (x=l, color='k', ls='dashed') 
pit . legend (loc=' best' ) 



pit . show ( ) 
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Ulesanne 73. d) 



Hinnake viga Taylori valemi kaudu saadud ligikaudses valemis: cosh a; w 1 

In [9] : def yl73d(x) : 

return np.cosh(x) 

def yl73d2 (x) : 

return 1 + x**2/2 

xl = get_x(0, 0.2) 
x2 = get_x ( 0 , 3 ) 

fig() 



kui \x\ < 0.2. 



pit . figure ( 1 ) 

pit . subplot (221) 

pit. plot (xl, yl73d(xl), 

pit. plot (xl, yl73d2(xl) 

pit . legend (loc=' best' ) 



r', label=' $\cosh{ (x) }$' ) 

'g', label=r'$l + \f rac { x"2 } { 2 } $ ' ) 



pit . subplot (222) 

pit. plot (x2, yl73d(x2), ' r' , label = ' $ \cosh { (x) } $ ' ) 

pit. plot (x2, yl73d2(x2), 'g', label=r' $1 + \f rac { x"2 } { 2 } $ ' ) 

pit . axvline (x=-0 . 2 , color='k', ls='dashed') 

pit . axvline (x=0 . 2 , color=' k' , ls='dashed') 

pit . legend (loc=' best' ) 



pit . show ( ) 




Ulesanne 73. e) 



Hinnake viga Taylori valemi kaudu saadud ligikaudses valemis: e x fal + x+ ^ + '^r-, kui |x| < 1. 



In [10] : def yl73e (x) : 

return e**x 

def yl73e2 (x) : 

return 1 + x + x**2/2 +x**3/6 

xl = get_x ( 0 , 1 ) 
x2 = get_x(0, 4) 

fig() 

pit . figure ( 1 ) 
pit . subplot (221) 

pit. plot (xl, yl73e(xl), ' r' , label=' $e A x$ ' ) 

plt.plot(xl, yl73e2(xl), ' g' , label=r ' $ l+x+\ f rac { x A 2 } { 2 } +\ f rac { x A 3 } { 6 } $ ' ) 
pit . legend (loc=' best' ) 

pit . subplot (222) 

pit. plot (x2, yl73e(x2), ' r' , label=' $e A x$ ' ) 

pit. plot (x2, yl73e2(x2), 'g', label = r ' $ l+x+\ f rac { x"2 } { 2 } +\ f rac { x A 3 } { 6 } $ ' ) 
pit . axvline (x=-l , color='k', ls=' dashed' ) 
pit . axvline (x=l , color='k', ls='dashed') 
pit . legend (loc=' best' ) 



pit . show ( ) 



